N87-22745 


0RK3SWK- WO- * 

OF POOR QOALJIi MODELING of controlled flexible structures WITH IMTULSIVE LOADS 


M. Zak 

Applied Technologies Section, Jet Propulsion Laboratory 
California Institute of Technology. Pasadena, CA, USA 


Abstract . The characteri stlc wave approach is developed as an alternative to 
modal methods which may lead to significant errors in the presence of 
impulsive or concentrated loads. The method Is applied to periodic 
structures. Some special phenomena like cumulation effects and transitions to 
ergodlclty are analyzed. 


INTRODUCTION 

Controlled large space structures, which will 
likely be composed of networks of long slender 
members, are subjected to disturbances (coming 
from the actuators) with a ralatively small 
contact zone and short time Interval. From the 
mathematical viewpoint such disturbances are 
characterized by discontinuities which can be 
considered as a very high frequency. In 
truncation techniques which are used in modal 
analysis tha contribution of high frequenclas Is 
lost, and thcrefora, the Impulsive concentrated 
loads are supposed to be treated by some other 
methods . 

Since any discontinuity propagates with the 
characteristic speed, it Is reasonable to turn to 
the characteristic wave approach In treating the 
Impulsive loads. The advantage of this approach 
Is In the fact that characteristic speeds depend 
only on the coefficients at the highest (second 
order) derivatives In the governing equation of 
structural members which significantly simplifies 
the analysis of characteristic waves. 

Thus, it appears that the application of the 
characteristic wave approach Is the most 
beneficial In the domains where spectral methods 
fall. That Is why It can be used as a supplement 
to modal methods for linear analysis of con- 
trolled structures when loads can be decomposed 
in to "smooth" and impulsiva components. 

In this article, some aspects of characteristic 
wave propagation, reflection and transmission in 
structures with one -dimensional structural 
■embers as well as possible engineering tools for 
their analysis are discussed. 


PROPAGATION OF IMPULSIVE LOADS IN 
ONE- DIMENSIONAL STRUCTURAL MEMBERS 


Here Cg, Cj, Cy, C M , and 0$ are the 
characteristic apeeds for longitudinal, 
torsional, shear, bending and transverse string 
waves, respectively. E is the Young modulus, G 
la the shear modulus, A is the cross-sectional 
area, A s la an effective shear area of a 
Timoshenko beam. El la the cross-sectional area 
moment of Inertia, pi Is the rotatory inertia per 
unit length, p la the mass per init length, T Is 
the string tension, U. Flugge, 1962. 

For homogenious structural members all the 
characteristic speeds (2) are constant, and 
consequently, the width Af as well as the 
duration At of the Impulse will be constant too. 
However, the Initial configuration of the impulse 
will be preserved only for the simple wave 
equation without damping. 

In all other cases due to the dispersion 
phenomenon this configuration, strictly speaking, 
will not be preserved. Nevertheless, the 
dispersion can be ignored If the conditions (1) 
are satisfied. For further convenience we will 
Introduce an equivalent rectangular Impulse of 
the same length and energy. For such a 
rectangular impulse, all the waves listed in (2) 
are decoupled even if they propagate 
simultaneously In a structural member, and this 
Is the qost Important advantage of the 
characteristic wave approach to propagation of 
Impulsive loads. 

As follows from the energy conservation law the 
height of the rectangular impulse expressed In 
terms of velocity, strain or stress will be 
constant If there is in material damping. If 
material damping Is proportional to the velocity 
being characterized by the damping coefficient f 
then the height h of the rectangular impulse will 
exponentially decrease: 

h 2 - h 2 0 (3) 


NON -HOMOGENEITY EFFECTS 


We will start with a one -dimensional structural 
mtmber subjected to a concentrated or impulsive 
load assuming thst 

At << L, and At << t 

C (1) 

in which L is the length of the structural 
member, Ai la the width of the contact zone of 
the Impulse, At is the duration of the 
concentrated load, and C is the characteristic 
speed of wave propagation, while 


The situation becomes more complicated even for a 
rectangular impulse if the speed of propagation 
Is not constant. This effect can be caused by 
non-llnear material properties (if the impulse is 
large enough to generate finite strains) or by 
non-homogenlous properties of the structural 
member. 

In the first esse the speed becomes non- 
characteristic since it depends on the magnitude 
of the transmitted parameters, i.e., for a 
rectangular impulse: 
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c - c(h) (4) 

This dependence may lead to a qualitatively new 
effects such as shock wave formation. 
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In the second esse the speed remains 
characteristic, but it depends on the spnee 
coordinate : 

c - c(x) (5) 

Although the governing equations for wave 
propagations remain linear (but with variable 
coefficients) the dependence (5) may lead to some 
surprising effects. In order to describe them, 
ve will start vlth the energy balance. The 
energy E of a propagating impulse consist of 
potential and kinetic components: 

E - i ( vf> dx <6) 
X 1 

in which xi and x 2 are the space coordinates of 
the trailing and leading fronts of the 
propagating wave, (v] and [<] are the impulsive 
velocity and strain, respectively, while 

it) - -c[v) (7) 


as follows from the klnematlcal condition at the 
front of a discontinuity, Miklovitz, (1984). 
Hence, for a rectangular Impulse: 

E - Hv) 2 (x 2 • X!) (8) 

and instead of Eq. (3) now one obtains: 

h 2 (x 2 - x i> - h 2 Aie“0 (9) 

The first effect which can be found from F.q. (9) 
Is associated with the specific energy 
cumulation, i.e., with the unbounded growth of h 
due to shrinking width ( x 2 - x^) of the 
propagating impulse. This effect was first 
described and explained by M. Zak, 1983. 


The second effect Is associated with a trapping 
of a propagating impulse within a localized area 
of a structural member. Similar effect of normal 
mode localization was predicted by C.H. Hodges, 
1962, in connection vlth a system of coupled 
linear oscillators with damping. Thus, the 
trapping effect of a propagating Impulse can be 
considered as a "continuum version" of the normal 
mode localization. A mathematical treatment of 
this effect is presented below. 


Consider a function (5) in the following form: 


C - 


C G at x < x^ and x > x^ 

7 C 0 at x^> x > x^, 0 < y < 1 


( 10 ) 


Such a discontinuity of the characteristic speed 
C within a small segment (*** - x^) < A f can be 
caused by some structural Irregularities (such as 
material inclusions, joints, etc.) 


As follows from Eq. (10) 


x 2 *xi 


Al at xj - x ^- Af , x 2 " x * 

Af-C 0 (l-7)t at x 1 <x^<x 2 <%* 

Al ' (x *** x *l x l <x ** x 2 >x ** (ID 

Af* (x #+ -x^)+C 0 (l- r)t at \'M <X **-X2 >x t * 
Af at xi , x 2 > x** 


Substituting OD into (9) one finds: 
hV ;t »t x t «x s^x^x. 

hV U /(l-Ct) at x <x„<x <x.„ 
o 1 

h : . { h|e' U /(l-» x l <x *l <x 2 <x ** 


hV'Vo-V^t) at x„<x <x x,> 
o * 


„2 - Ct 
h,e 


y at x., X 2 > x. 


in which 


c ( i-O 


^ = At “ AO ' 


4D 


2 2 -tt 2 2 t< x* t 2 ) 

• h l * V • h 2 ■ V 


2 t 2 '"Wh’ 


h‘j • h Q e 


. x . L - r-x , x 

h ■ r 2 * c • 3 * i 


Simple analysis of Eq . (1?) shows 

function h(x) has maximum at xj < < 

if 


< < 1 < & 


This maximum is: 

h ■ n A = e 
max 0 t. 


The maximum of h max as a function of r / 

at 


t 

l 


| (1 ♦ /5) 


and therefore: 

h2 max.max * 8.15 h 2 0 
Substituting (16) into (14) one obtains 
X < 0.382 


But, as follows from (13): 


.2 .2 
h, * h e 
3 o 


- | (l4 Y **)t 


Hence, the trapping effect will be the 
if 


X - 0.382 


02 ) 


03) 

thnt the 

x 2 < x ** 

(14) 

(15) 
f will be 

(16) 

07) 

OB) 

(19) 
strongest 

(20) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Indeed, In this esse the function (12) has the 
sharpest shape since Eq. (16) provides the 
largest maximum (17), while Eq . (20) leads to the 
highest degree of dissipation after this maximum. 
Thus, the conditions (16) and (20) can be used as 
the key for structural implementation of the 
trapping effect. 


Governing Equation* 


As shown by Zak, M. 1965, the governing equations 
for a propagation, reflection and transmission of 
an initial impulse In a two -member structure with 
Isolated ends at x - 1 , and x - 3, and a Joint at 
x - 2 form a system of difference equations: 

h k4l * Ah k < 21) 

where 


o run 0 0 

-fl2f2 z 2 0 0 -ri2f2< 1 - z 2 

-f23f2< 1 ' z 2> 0 0 "f 23^2 Z 2 

0 0 -r23f3 0 



( 22 ) 


in which h }2 Is the wave at x • 2 coming from x - 
1 , etc. f 12 and £33 are the damping coefficients 
in the corresponding structural members, f*, 

{■3 are the damping coefficients at x - 1 , 2 , 3 , 
respectively, ^2 * s th * reflection coefficient at 
X - 2, while t x . t 3 - 1 . 

The same procedure can be applied to multi -member 
structural systems, while the matrix A will 
attain new submatrices corresponding to 
additional structural hers (like a stiffness 
man ix in finite element methods) . 

Thus, ary structure with n identical structural 
members subjected to impulsive loads is governed 
by the matrix difference equation ( 21 ) while the 
order of this system is 2 n. 

AmlYSls oflSolutlons 

In the case of n-member structure the solution to 
the governing equation (21) where the matrix A Is 
of the order 2 n can be written in the following 
form: 


h k - A k h e , k- 0.1.2.... (23) 

All the qualitative properties of this solution 
are defined by the eigenvalues of the matrix A, 
by the roots of the characteristic 
polynomial : 

|A • Al| - 0 (24) 

For Instance, the solution is stable if 

l«l <1.1- 1.2. ... 2n (25) 


By applying a linear fractional transformation 


A « 


X-l 


(26) 


to Eq. (25) one reduces the stability analysis to 
conventional methods. 


Transition to Ereodicltv 

So far Lhe only structures with idem Icq] members 
(characterized by the dlmentlonless time delay I) 
were considered. It was demonstrated that there 
exists a formal analogy between the matrix 
techniques for treating these structures under 
Impulsive loads and for conventional modal 
anilysls (although the matrices have different 
physical nature). However, in reality the 
Identicalness of structural member* even in 
periodic structures 1* an exception rather than a 
rule. Indeed, different time delays for 

different structural members can be caused not 
only by different lengths, but also by different 
characteristic speeds. In turn, different 
characteristic speeds may i>ccur if joints convert 
one type of deformation into another (see Eq. 
(2)). Another source of different time delays is 
associated with external forces If they applied 
not to joints. In this case, the points of their 
application must be considered as additional 
Joints, but without reflection or damping, and 
this will lead to additional structural members 
with different time delays. As will be shown 
below, different time delays lead to new 
qualitative effects which do not occur in modal 
methods. 

For simplicity, we will start with the tvo*membc: 
structure and assume the following 
(dimensionless) lengths: 

AB - 4, BC - 6 (27) 

Then, the characteristic equation for this t 
member system obviously has the order 24 (which 
is the least common multiple of 2x4 and 2x6) 

But In the case 

AB - 1, BC - (28) 

the ratio of the time delays is irrational, and 
therefore, for successive rational approximations 
of 77 the order of the governing difference 
equation tends to infinity as: 

4, 14, 282, 1414, ... etc. (29) 

Now it Is easy to deduce, that if in an n-member 
structure the time delays rj, r 2t ... r n .. . c 
commensurate, then the order of the governing 
difference equation will be finite and equal to 
the least common multiple of 2 ri, 2 r 2 , 2 * n . 

If at least two time delays are not coramersurate , 
this order will tend to infinity. Obviously, 
this effect does not have an analogy in modal 
approach where the order of the governing 
differential equation depends only on the number 
of modes (or finite elements) considered. 


In order to clarify the physical meaning of such 
a phenomenon let us start with the following 
question: during what time Interval T an initial 

Impulse will return to its original location in 
"one piece"? Simple geometrical consideration 
show that 

T - 2 if AB - 1. BC - 1 

T - 24 if AB - l, BC - 6 (30) 

T - • If AB - 1, BC - y? 


In other words, this (dimensionless) interval is 
equal to the order of the governing difference 
equation . 
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Thus, if at least two tine delays in an n-nember 
structure are not connensurate, the system will 
never return to its initial position, i.e, t the 
notion will lose its periodicity. In classical 
mechanics such systems are known as ergodlc 
systems. For infinite number of tines they pass 
through every. state of notion (which Is 
consistent with constraints) spending equal tine 
Intervals near ech state. 

As however, the rational numbers are a set of 
measure aero, practically every notion of this 
type sooner or later become ergodlc. 

Nevertheless, in engineering applications there 
always can be found such a characteristic time 
interval within which the notion is approximately 
periodic, while the transition to ergodicity can 
be Ignored due to damping. 

It is worth emphasizing that the transition to 
ergodicity is not •inevitable" if one takes into 
account non* linear properties of real structures. 
Non-linearities nay provide some mechanisms (such 
as dynamical synchronization effects) which 
depress the disorder and lead to periodical 
notion. In this connection, it is relevant to 
mention the experiment with coupled chain of 
harmonis oscillators performed by Fermi, Pasta 
and Ulam. Instead of ergodicity which was 
expected they found periodic oscillations. 
However, if dynamical synchronization effects do 
not depress ergodicity end if the characteristic 
time during which the notion can be approximated 
as periodic is too short one has to apply methods 
of statistical mechanics. 
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